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Abstract: We investigate a supersymmetric extension of the Minimal Supersymmetric 
Standard Model (MSSM), called the TNMSSM, containing a SU{2) Higgs triplet (T) of 
Y = 0 hypercharge and a singlet superfields ( S ) in the corresponding superpotential. The 
model can be viewed, equivalently, as an extension of the NMSSM with the addition of 
a T — S interaction and of an extra coupling of the triplet to the two Higgs doublets 
of the NMSSM. In this scenario the Higgs particle spectrum at tree-level gets additional 
mass contributions from the triplet and singlet scalar components respect to the MSSM, 
which are particularly enhanced at low tan f3. We calculate the one-loop Higgs masses 
for the neutral physical Higgs bosons by a Coleman-Weinberg effective potential approach. 
In particular, we investigate separately the impact of the radiative corrections due to the 
electroweak, gauge-gaugino-higgsino, fermion-sfermion and Higgs self-interactions to the 
Higgs masses. Due to the larger number of scalars and of triplet and singlet couplings, the 
Higgs corrections can be larger than the strong corrections. This reduces the amount of 
fine-tuning required to fit the recent Higgs data. Using the expressions of the beta-functions 
of the model, we show that the large triplet singlet coupling remains perturbative up to 
~ 10 8 ~ 10 GeV. The model is also characterized by a light pseudoscalar in the spectrum, 
which is a linear combination of the triplet, doublet and singlet CP-odd components. We 
discuss the production and decay signatures of the Higgs bosons in this model, including 
scenarios with hidden Higgses, which could be investigated at the LHC in the current run. 
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1 Introduction 

With the recent discovery of the Higgs boson at the Large Hadron Collider, the mechanism 
responsible for the breaking of the electroweak symmetry has finally been uncovered and it 
has been shown to involve at least one scalar field along the lines of the Standard Model 
(SM) description. This discovery has removed, at least in part, previous doubts about the 
real existence of a scalar with Higgs-like properties in our Universe. Both the CMS [1, 2] 
and the ATLAS [3] experimental collaborations have confirmed the discovery of a Higgs 
boson, by an analysis of the 77 , ZZ*, and WW* decay channels of the Higgs particle - as 
predicted by the Standard Model (SM) - at a confidence level of more than 5a, except for the 
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WW* decay rate, which has been recorded with a 4.7 a accuracy by CMS [2], The fermionic 
decay modes, instead, have still to reach the 5a accuracy, and show some disagreement in 
the results elaborated by the two experimental collaborations. Clearly, the disagreement of 
the experimental results with the predictions from the SM opens the possibility of further 
investigation of the Higgs sector. 

For such reasons, it is widely believed that the SM is not a complete theory, being not 
able, for instance, to account for the neutrino masses, but also for being affected, in the 
scalar sector, by the gauge hierarchy problem [4], The widespread interest in the study of a 
possible supersymmetric extension of this model has always being motivated with the goal 
of finding a natural and elegant solution to this problem. In fact, supersymmetry protects 
the Higgs mass from the undesired quadratic divergences introduced by the radiative cor¬ 
rections in the scalar sector of the SM, and it does so by the inclusion of superpartners. 

In the minimal supersymmetric extension of the SM (MSSM) the conditions of analyticity 
of the superpotential and of absence of the gauge anomalies require a minimal extensions 
of the scalar sector with two Higgs superfields, in the forms of SU( 2) doublets carrying 
opposite hypercharges (Y). Supersymmetric extensions are, in general, characterized by 
a large set of additional parameters which render their phenomenological study quite in¬ 
volved. For this reason, in the near past, the interest has turned towards models, such as 
the constrained minimal supersymmetric extension (cMSSM/mSUGRA), with only 5 new 
parameters, generated at a large supergravity scale, quite close to the Planck scale [5]. 

Unlike the SM case, in the MSSM the tree-level mass of the lightest Higgs (hi) m 
is not a free parameter but it is constrained to lay below the mass of the Z gauge boson, 
mz (m/i! < rriz ) ■ This constraint has been in tension with the results of the experimental 
searches at LEP-2 which have failed to detect any CP-even Higgs below mz and which 
had established a lower bound of 114.5 GeV for the SM Higgs boson [6]. With the recent 
discovery of a CP-even Higgs boson around 125 GeV [1-3] the resolution of this conflict is, 
therefore, mandatory. 

Even in this unfavourable situation, supersymmetric scenarios remain popular due their 
naturalness and for having a dark matter candidate in theories with a conserved i?-parity 
symmetry [7]. To avoid the conflict between the MSSM prediction for the Higgs and the 
LHC results, one needs to consider the effect of the radiative corrections which could lift 
the bound on the Higgs mass in this model. It has been shown - and it is now well known - 
that in the case of the MSSM the significant radiative corrections come from the stop-top 
corrections, specially at low tan /?, due to large Yukawa couplings and to the presence of 
colour charges. This has triggered analysis envisioning scenarios with a heavy stop, which 
require a very high supersymmetric (SUSY) scale for the most constrained supersymmetric 
models like mSUGRA/cMSSM, AMSB, etc [8]. In the case of the phenomenological MSSM 
(pMSSM) there are two possibilities: a very large third generation SUSY mass scale and/or 
a large splitting between the stop mass eigenstates [9]. The second case leads to large soft 
trilinear couplings > 2 TeV [9], which brings back the fine-tuning problem in a different 
way. 

A possible way to address the fine tuning problem is to consider an extended Higgs sector. 
In this respect, there are some choices which could resolve it, based on the inclusion of 
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one singlet [10] and of one or more triplet superfields of appropriate hypercharges [11]. In 
particular, the addition of a Y = 0 hypercharge superfield gives large tree-level as well as 
one-loop corrections to the Higgs masses, and relaxes the fine tuning problem of the MSSM 
by requiring a lower SUSY mass scale [12, 13]. There are some special features of these 
extensions which are particularly interesting and carry specific signatures. For instance, the 
addition of a (Y = 0, ±2 hypercharge) - (SU( 2) triplet) Higgs sector induces H± — W T — Z 
couplings mediated by the non-zero vacuum expectation value (vev) of the Higgs triplet, due 
to the breaking of the custodial symmetry [14, 15]. Other original features of the Y = ±2 
hypercharge triplets are the presence of doubly charged Higgs in the spectrum [16] . There 
are also other significant constraints which are typical of these scenarios, and which may 
help in the experimental analysis. In the supersymmetric Higgs triplet extension, the vev 
of the triplet vt is highly constraints by the p parameter [17], which leads to vt < 5 GeV 
in the case of Y = 0 triplets. In the same case, this value of vt can account for the value of 
the mixing parameter ^ of the 2-Higgs doublets (or /x-term), which remains small in the 
various possible scenarios. Another dynamical way to generate a po term is by adding a 
SM gauge singlet superfield to the spectrum [10], as in the NMSSM. Thus a triplet-singlet 
extended supersymmetric SM built on the superpotential of the MSSM, can address both 
the fine tuning issue and resolve, at the same time, the problem of the /x-term of the two 
Higgs doublets [18-20]. 

In this work we are going to investigate these extensions, presenting results on the spectrum 
of these models. We will also see that the addition of a discrete symmetry in this model 
removes the mass terms from the super potential and its continuum limit generates a Nambu- 
Goldstone pseudoscalar particle in the spectrum, characterising some of its most significant 
features. A more thorough analysis of this specific aspect will be presented elsewhere. The 
goal of our study is to investigate the allowed region of their parameter space in view of 
the experimental constraints emerging from the recent experimental results at the LHC. 
Our work is organized as follows. We start by introducing a scale invariant superpotential 
which is MSSM-like, but with the inclusion of a Y = 0 Higgs triplet and of an extra SM 
gauge singlet superfield. In section 2 we detail the model. In section 3 we investigate 
both the possibility of generating a tree-level Higgs mass around 125 GeV due to the extra 
contributions from the triplet and the singlet, and the possibility of hidden Higgs bosons. 
Section 4 describes the strong and weak sectors of the model, followed by a discussion 
of the strong and weak contributions to the radiative corrections of the Higgs masses in 
section 5. We address in section 6 the perturbative running of the couplings as we vary 
the scale of the theory, and in section 7 we examine the issue of fine-tuning in this model. 
A study of the light pseudoscalar state is investigated in section 8. Finally, in section 9 
we take into account the constraints from the LHC Higgs data and LEP data to find out 
the phenomenological parameter space. In section 10 the new possible production and 
decay channels are given along with the signatures that could be tested at the LHC and 
at the future colliders. In section 11 we present our conclusions, summarizing our results 
and elaborate on possible extensions of this work. We have left to an appendix both the 
expressions of the renormalization group equations of the dimensionless couplings and the 
tree-level vertices of the Higgs sector. 
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2 The Model 


We consider a scale invariant super potential Wtnmssm with an extended Higgs sector 
containing a. Y = 0 SU(2) triplet T and a SM gauge singlet S (see [18, 19]) on top of the 
superpotential of the MSSM. We recall that the inclusion of the singlet superfield on the 
superpotential of the MSSM realizes the NMSSM superpotential. We prefer to separate the 
complete superpotential of the model into a MSSM part, 


Wmssm = UtUHu-Q - UbDHd-Q ~ UrEHd-L , (2.1) 


where denotes a contraction with the Levi-Civita symbol , with e 12 = +1, and combine 
the singlet superfield (5) and the triplet contributions into a second superpotential 


W TS = A T H d • TH U + X s SH d -H U +^S 3 + X ts STt[T 2 } 


with 


Wtnmssm = Wmssm + Wts- 
The triplet and doublet superfields are given by 


( 2 . 2 ) 

(2.3) 


t = 





(2.4) 


Here T° is a complex neutral superfield, while and are the charged Higgs superfields. 
Note that ( T , 1 - )* ^ ■ Only the MSSM Higgs doublets couple to the fermion multiplet 

via Yukawa coupling as in Eq. (2.1), while the singlet and the triplet superfields generate 
the supersymmetric term after their neutral parts acquire vevs, as shown in Eq. (2.2). 

In any scale invariant supersymmetric theory with a cubic superpotential, the complete 
Lagrangian with the soft SUSY breaking terms has an accidental Z 3 symmetry, the invari¬ 
ance after the multiplication of all the components of the chiral superfield by the phase 
e 2m / 3 . Such terms are given by 


Woft — m 2 H u \H u \ 2 + m Ti d \Hd\~ + mils'! 2 + rnr\T\ 2 + mg|Q|'+ m^-|?7| 2 + mJj\D\~ 
+(A s SH d .H u + + A T H d .T.H u + A TS STr(T 2 ) 

+ AuUHu-Q+ AdDHd.Q + h.c), (2-5) 


while the D-terms are given by 

Vb = \ tij < Ai) 2 - (2■-6) 

k 


In this article we assume that all the coefficients involved in the Higgs sector are real 
in order to preserve CP invariance. The breaking of the SU(2)l x £/(l)y electroweak 
symmetry is obtained by giving real vevs to the neutral components of the Higgs fields 


^ IT'J u 

u V2’ 


<H° d 


>=^=, ,<S>=^= 

Vi Vi 


< T °> = v I’ (2,7) 
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which give mass to the W ^ and Z bosons 

™w = \aHv 2 + 44), m% = |(4 + 9y)v 2 , v 2 = (4 + 4)- 

and also generate the p£, = ^vs + ^~vt term. 

The non-zero triplet contribution to the W± mass leads to a deviation of the tree-level 
expression of the p parameter 


P ~ 1 + 4 777 • 


(2.9) 


Thus the triplet vev is strongly constrained by the global fit on the measurement of the p 
parameter [17] 

p= 1.00041™, (2-10) 

which restricts its value to Vt < 5 GeV. In our numerical analysis we have chosen vt = 3 
GeV. 

3 Tree-level Higgs masses 

To determine the tree-level mass spectrum, we first consider the tree-level minimisation 
conditions, 

d^V\ vev = 0; V = V D + V F + V soft , < >= -*|, T,; = H°,H° d ,S,T 0 , (3.1) 

where we have defined the vacuum parameterizations of the fields in the Higgs sector as 

H° u = flS = -L(flS )r +iflS ii ), S = ±(S r +iSi), T° = -^(T°+zT°). 

(3.2) 

from which the soft-breaking masses are derived in the form 


m H = 


Vd 
2 v u 


m H d = 


i^/2Asvs — vt (yAr + VZvsXt^ts'J + As (nvg + n^Ars)) 

— 2 ^ ~ 4) + 2^ (4 + 4) + V2\svs\tvt^J 

+ g (4 _ 4)(4 + 4), ( 3 - 3 ) 

(\/iT4sns - n T + As (ku| + 4Ats)) 

_ 2 (a| (4 + 4) + 2^t (4 + 4) - V^XsvsXtvt'J 
+ o (4 _ 4) (4 + 4), ( 3 - 4 ) 


m s = 


1 


242 ns 

-4 K ns 


(nr (At (As (4 + 4) — ^ v d v u^Ts) - 2 Atsvt) + 2AsVdV u ) 


+ KVdV u \s — —As (n^ + 4) — K 2 n| — kvtXts ~ 2 v t^tSi (3-5) 
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m T = ~^~ (V%vsXt (As (v d + v^) - 2v d v u X ts) ~ 2 AtvaVu^J - \[2Atsvs 

+ A ts {v d v u \ s — vg (k + 2X ts)) ~ ^Aj. (t^ + — Vj-X ts- (3-6) 

It can be shown that the second derivative of the potential with respect to the fields 
satisfy the tree-level stability constraints. The neutral CP-even mass matrix in this case is 
4-by-4, since the mixing terms involve the two 517(2) Higgs doublets, the scalar singlet 5 
and the neutral component of the Higgs triplet. After electroweak symmetry breaking, the 
neutral Goldstone gives mass to the Z boson and the charged Goldstone bosons give mass 
to the boson. Being the Lagrangean CP-symmetric, we are left with four CP-even, 
three CP-odd and three charged Higgs bosons as shown below 

CP — even CP — odd charged 
hi,h2,h 3 ,h 4 01 , 02,03 hf,hf,hf. (3.7) 

The neutral Higgs bosons are combination of doublets, triplet and singlet, whereas the 
charged Higgses are a combination of doublets and triplet only. We will denote with m di 
the corresponding mass eigenvalues, assuming that one of them will coincide with the 125 
GeV Higgs (/ 1125 ) boson detected at the LHC. The scenarios that we consider do not assume 
that this is the lightest eigenvalue which is allowed in the spectrum of the theory. Both 
scenarios with lighter and heavier undetected Higgs states will be considered. In particular, 
we will refer to those in which one or more Higgses with a mass lower than 125 GeV is 
present, to hidden Higgs scenarios. 

At tree-level the maximum value of the lightest neutral Higgs has additional contributions 
from the triplet and the singlet sectors respectively. The numerical value of the upper 
bound on the lightest CP-even Higgs can be extracted from the relation 

ml, < m|(cos 2 2/3 + 9 ^ 9 sin 2 2/3 4 - 9 9 sin 2 2/3), tan/3 = —, (3.8) 

91 + 9y 9 l + 9y y d 

which is affected on its right-hand-side by two additional contributions from the triplet 
and singlet. These can raise the allowed tree-level Higgs mass. Both contributions are 
proportional to sin2/3, and thus they can be large for a low value of tan/3, as shown in 
Figure 1. The plots indicate that for higher values of A t,s a lightest tree-level Higgs boson 
mass of ~ 125 GeV can be easily achieved. For general parameters, the required quantum 
corrections needed in order to raise the mass bound are thus much smaller compared to the 
MSSM. In the case of the MSSM, as we have already mentioned, at tree-level < mz, 
and we need a correction > 35 GeV to match the experimental value of the discovered 
Higgs boson mass, which leads to a fine-tuning of the SUSY parameters. In fact, this 
requires that the allowed parameter space of the MSSM is characterized either by large 
SUSY masses or by large splittings among the mass eigenvalues. In fact, this requires that 
the allowed parameter space of the MSSM is characterized either by large SUSY masses 
or by large splittings among the mass eigenvalues. We have first investigated the tree-level 
mass spectrum for the Higgs bosons and analysed the prospect of a ~ 125 GeV Higgs boson 
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m hl 



Figure 1. Tree-level lightest CP-even Higgs mass maximum values with respect to tan/3 for (i) 
A t = 0 . 8 , A s = 0.1 (in red), (ii)AT = 0 . 1 , A s = 0.8 (in green ) and (iii) X T = 0 . 8 , A s = 0.8 (in 
blue). 


along with the hidden Higgs scenarios. We have looked for tree-level mass eigenvalues where 
at least one of them corresponds to the Higgs discovered at the LHC. For this purpose we 
have performed an initial scan of the parameter space 

< 1, |^| < 3, |u s | < 1 TeV, 1 < tan/3 < 10, (3-9) 

and searched for a CP-even Higgs boson around 100 — 150 GeV, assuming that at least one 
of the 4 eigenvalues will fall within the interval 123 GeV < rrih- < 127 GeV at one-loop. 

Figure 2(a) presents the mass correlations between rri} M and rrih 2 , where we have a CP- 
even neutral Higgs boson in the 100 < mh. t < 150 GeV range. The candidate Higgs boson 
around 125 GeV will be determined at one loop level by including positive and negative 
radiative corrections in the next section. The mass correlation plot at tree-level shows that 
there are solutions with very light h\, m./ ll < 100 GeV, which should be confronted with 
LEP data [6]. At LEP were conducted searches for the Higgs boson via the e + e~ —> Zh 
and e + e _ —> h\h 2 channels (in models with multiple Higgs bosons) and their fermionic 
decay modes (h —> bb,rf and Z —> t€). The higher centre of mass energy at LEP II (210 
GeV) allowed to set a lower bound of 114.5 on the SM-like Higgs boson and of 93 GeV 
for the MSSM-like Higgs boson in the maximal mixing scenario [ 6 ]. Interestingly, neither 
the triplet (in our case) nor the singlet type Higgs boson couple to Z or to leptons (see 
Eq. (2.2)), and as such they are not excluded by LEP data. 

We mark such points with > 90% triplet/singlet components, which can evade the LEP 
bounds, in green. In Figure 2(a) one can immediately realize that the model allows for some 
very light Higgs bosons (m/ il < 100 GeV). We expect that the possibility of such a hidden 
Higgs would be explored at the LHC with 14 TeV centre of mass energy, whereas the points 
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where h\ is mostly a doublet (> 90%) could be ruled out by the LEP data. The points 
with the mixed scenario for h\ (with doublet, triplet and singlet) are marked in blue. We 
remark that a triplet of non-zero hypercharge will not easily satisfy the constraints from 
LEP, due to its coupling to the Z boson. 

For the points with < 100 GeV which are mostly doublet (red ones) it is very 

hard to satisfy the LEP bounds [ 6 ]. This is because, being doublet like, such h\ would 
have been produced at LEP and decayed to the fermionic pairs, which have been searched 
extensively at LEP. On the other hand the singlet and triplet like points (green points) are 
very difficult to produce at LEP due to the non-coupling to Z boson, which was one of the 
dominant production channel. This is true for both e + e _ —> Zh\ and e + e _ —> h\a\. Such 
triplet and singlet like points will reduce the decay widths in charged lepton pair modes 
due to non-coupling with fermions. These make the green points more suitable candidate 
for the hidden Higgs bosons, both for the CP-even and CP-odd. However such parameter 
space would be highly constrained from the data of the discovered Higgs boson around 
125 GeV at the LHC. So far the discovery of the Higgs boson at the LHC has reached 
5 <t or more in the channels /1125 —> 77 , WW*, ZZ*. Effectively this could be satisfied by 
the candidate Higgs around 125 GeV which is mostly doublet like and its decay branching 
fractions should be within the uncertainties give by CMS and ATLAS experiments at the 
LHC. Such requirements rule out vast number of parameter points, including some the 
triplet and/or signet like hidden Higgs boson(s). In section 9 we consider such constrains 
coming from the Higgs data at LHC and the existing data from LEP. 

Figure 2(b) shows the mass correlation between I 13 and /14 for the the same region (3.9) 
of the parameter space. We see that although there are points characterized by a mass rrih 3 
lighter than 500 GeV, states with m / l4 < 500 GeV are less probable. Figure 3 shows the 
mass correlations of the CP-odd neutral Higgs bosons. Specifically, Figure 3(a) presents the 
analysis of the mass correlation between 07 and 02 - The plot shows that there exists the 
possibility of having a pseudo-scalar ai lighter than 100 GeV, accompanied by a CP-even 
~ 125 GeV Higgs boson. Note that a very light pseudoscalar Higgs in the MSSM gets strong 
bounds from LEP [ 6 ]. In this case, for a high tan/3, the pair production process e + e _ —>• hA, 
where A is the pseudoscalar of the MSSM, is the most useful one, providing limits in the 
vicinity of 93 GeV for 7774 [ 6 ]. In the TNMSSM instead, if the light pseudoscalar Higgs 
bosons are either of triplet or singlet type then they do not couple to the Z, which makes it 
easier for these states to satisfy the LEP bounds. For this purpose, points which are mostly- 
triplet or -singlet (90%) have been marked in green; points which are mostly-doublet (90%) 
in red, whereas the mixed points have been marked in blue as before. Certainly, mass 
eigenvalues labelled in green would be much more easily allowed by the LEP data, but 
they would also be able to evade the recent bounds from the LHC Htt decay mode for a 
pseudoscalar Higgs [21]. This occurs because neither the triplet nor the singlet Higgs boson 
couple to fermions (See Eq. (2.2)). Figure 3(b) presents the correlation between ao and 03 
where the same colour code applies for the structure of 02 - As one can easily realize from 
the figure, there are plenty of green coloured points which represent triplet/singlet type 02 
states, which can easily evade the recent bounds on pseudoscalar states derived at the LHC 
[21]. Figure 4 shows the correlation of the three charged Higgs bosons for the region in 
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Figure 2. Tree-level CP-even Higgs mass correlations (a) vs m h2 and (b) vs to/j 4 , where 
we have a candidate ~ 125 GeV Higgs boson. The colors refers to the character of the hi mass 
eigenstate, describing the weights of the doublet, singlet and triplet contributions in their linear 
combinations. Red points are > 90% doublets-like, the green points are either > 90% triplet-like 
or singlet-like and blue points are mixtures of doublet and triplet/singlet components. The linear 
combinations corresponding to green points are chosen to satisfy the constraints from LEP onto Z 
and lepton final states. 



Figure 3. Tree-level CP-odd Higgs mass correlations (a) m ai vs m a2 and (b) m a2 vs m a3 , where 
we have a candidate ~ 125 GeV Higgs boson. The red points are > 90% doublets-like and the green 
points are > 90% triplet-like. The blue points are mixtures of doublet and triplet components for 
ai in (a) and for a 2 in (b) respectively. 


parameter space where we can have a ~ 125 GeV Higgs candidate. Figure 4(a) shows that 
there are allowed points for a charged Higgs of light mass (jn h ± < 200 GeV) correlated with 
a heavier charged Higgs hf. Only Higgses of doublet and triplet type can contribute to 
the charged Higgs sector. We have checked the structure of the lightest charged Higgs hf 

















Figure 4. Tree-level charged Higgs mass correlations (a) m h ± vs rn h ± and (b) m h ± vs m h ±, 
where we have a candidate ~ 125 GeV Higgs boson. The red points are > 90% doublets-like and 
the green points are > 90% triplet- or singlet-like. The blue points are mixture of doublet and 
triplets/singlets, for hf in (a) and for hf in (b) respectively. 


in Figure 4(a), where the red points correspond to > 90% doublet, while the green points 
correspond to > 90% triplet and the blue points to doublet-triplet mixed states. Charged 
Higgs bosons which are mostly triplet-like in their content (the green points) do not couple 
to the fermions (see Eq. (2.2)), and thus can easily evade the bounds on the light charged 
Higgs derived at the LHC from the H± —>• tv decay channel [22], This kind of triplet 
charged Higgs boson would also be hard to produce from the conventional decay of the top 
quark and the new production modes as well as the decay modes will open up due to the 
new vertex hf — Z — W T [15]. Thus vector boson fusion (VBF) with the production of a 
single charged Higgs is a possibility due to a non-zero hf — Z — W T vertex [15]. Apart from 
the hf -4 ZW ^ channels, the hf -4 a\{h\)W± channels are also allowed, for very light 
neutral Higgs bosons ( a\/h \). Figure 4(b) presents the correlation between m,± and m,±. 

i ... 2 3 

We have used for hf the same colour conventions as in the previous plots. We see that 
there are only few triplet type hf (green points), most of the allowed mass points being 
doublet-triplet mixed states (blue points). 

4 Strong and weak sectors 

The TNMSSM scenario has an additional triplet which is colour singlet and electroweak 
charged and a singlet superfields (see Eq. (2.2)) not charged under SU (3) c X SU (2 )l x U (l)y. 
Therefore, the strong sector of the model is the same of the MSSM, but supersymmetric 
F-terms affect the fermion mass matrices, and contribute to the off-diagonal terms. It 
generates additional terms in the stop mass matrix from the triplet and singlet vevs, which 
will be shown below. These terms are proportional to A tVt and Agug respectively, and 
allow to generate an effective /i£>-ternr in the model. The triplet contribution is of course 
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restricted, due to the bounds coming from the p parameter [17]. Thus, a large effective 
term can be spontaneously generated by the vev of the singlet, vs- 

Figure 5 shows the mass splitting between the t ,2 and t\ stops versus \g, for several vs 
choices and with At = 0. Large mass splittings can be generated without a large parameter 
At, by a suitably large vs, which is a common choice if the singlet is gauged respect to an 
extra U(l)' [23], due the mass bounds for the additional gauge boson Z' [24], The mass 
matrices for the stop and the sbottom are given by 



Figure 5. The mass splitting between the stop mass eigen states (£ 2 , 1 ) vs As for A t = 0 with 
vs = 500,1000,2000 GeV respectively. 


1 m b + m Q, + A (Sy + 3 9l ) ( v l ~ v i) 

V 7S Am + - »sA s ) 


m b + m % + h( v l- v d ) 9y 


\ 

(4.2) 


In the electroweak sector the neutralino (x^ =1 6 ) and chargino (x^i 2 3 ) sec t° r are 
enhanced due to the extra Higgs fields in the super potential given in (2.2). The neutralino 
sector is now composed of B, W 3 , H U1 H^, To, S. The corresponding mass matrix is thus 
now 6-by-6 and given by 
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/ M, 

0 -\gyv d 

\gyv u 

0 

0 

0 

m 2 \gLVd 

-\gLv u 

0 

0 

-\gyvd 

\9LVd 0 

\vt\t - ~^ v s^s 

\v u X T 

~72 VuXs 

\gYv u 

-\gLv u \vt^t - tjUsAs 

0 

\vdXr 

2 VdXs 

0 

0 \v u X T 

\v d \r 

V^vsXts V2vtXts 


\ 0 0 ~ 72 VuXs ~72 VdXs V^vtXts V2kv s . / 


(4.3) 

The triplino (To) and the singlino (S') masses and mixings are spontaneously generated by 
the corresponding vevs. The triplino and singlino are potential dark matter candidates and 
have an interesting phenomenology as they do not couple directly to the fermion superfields. 
The doublet-triplet (singlet) mixing is very crucial in determining the rare decay rates as 
well as the dark matter relic densities. 

Unlike the neutralino sector, the singlet superfield does not contribute to the chargino 
mass matrix, and hence the MSSM chargino mass matrix is extended by the triplets only. 
The chargino mass matrix in the basis of W + , H^, ^2 , Tj“) takes the form 


M, ± 


i M 2 ^giv u ~9 lvt 

^75 9LVd ^ ^vs\s + \vt\t 


(4.4) 


\ 9LVT - ^ v d^T \/2 VsXtS ) 

The chargino decays also have an interesting phenomenology due to the presence of a 
doublet-triplet mixing. 


5 Higgs masses at one-loop 


To study the effect of the radiative correction to the Higgs masses, we calculate the one-loop 
Higgs mass for the neutral Higgs bosons via the Coleman-Weinberg effective potential given 
in Eq. (5.1) 


kcw 


647T 2 


STr 





(5.1) 


where At 2 are the field-dependent mass matrices, /J, r is the renormalization scale, and the 
supertrace includes a factor of ( — \) 2 J (2J+\) for each particle of spin J in the loop. We have 
omitted additional charge and colour factors which should be appropriately included. The 
corresponding one-loop contribution to the neutral Higgs mass matrix is given by Eq. (5.2) 


(A MDa = 


<9 2 Vcw(^) 

d^id^j 


$ij dVcw($) 

<$i> d§i 


E l dm 2 k dm\ m\ 

k 32vr 2 <9<f>i d&j H /.1 2 

_ \ ^ 1 2 hi I 

32vr 2 k (4 >i) d®i \ 


+ E 


32vr 2 


m 2 k 



vev 


d2m l A ™£_ -A 

V 9r J vev 

®i,j = H° a rl H^ r , S r ,T® . (5.2) 
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Here, m | is the set of eigenvalues of the field-dependent mass matrices given in the 
equation above, and we remind that the real components of the neutral Higgs fields are 
defined as 


K = H d = T 2 W.r + «$,(). 

i(S r + iSi), = + (5.3) 

For simplicity we drop the supertrace expressions in Eq. (5.2), but for each particle 
the supertrace coefficient should be taken into account. 

Having characterized the entire sector of the TNMSSM, we gear up for the numerical 
evaluation of the one-loop neutral Higgs masses in the model. We have already seen in 
Eq. 3.8 that for low tan (3 the contribution of the radiative corrections required in order 
to reach the ~ 125 GeV Higgs mass, overcoming the tree-level bound in (3.7), is reduced. 
This is due to the additional Higgs and higgsinos running in the loops. In our analysis we 
have chosen the following subregion of the parameter space 

< E |ac| < 3, |n s | < ITeV, 1 < tan/3 < 10, 

\At,s,ts,u,d\ < 500, 14,1 < 1500, m Q 3 ,u3,d 3 < 1000, (5-4) 

65 < |Mi, 2 1 < 1000, 

that we have used in the computation of the Higgs boson mass. In this scan, we have 
included the radiative corrections to the mass eigenvalues at one-loop order of the neutral 
sector and retained only those sets of eigenvalues which contain one 125 GeV CP-even 
Higgs. We have selected the range 65 < |Mi, 2 | < 1000 in order to avoid the constraints on 
the Higgs invisible decay and use = 500 GeV for the numerical calculation. 

Figure 6 shows the radiative corrections to rri} M as A = m ^ loop — plotted 

against (a) A t, (b) A s and (c) k respectively. The red points show the corrections to 
m/ ll from the strong sector, due to the contributions generated by top-stop and bottom- 
sbottom running in the loops. The blue points include the corrections from the weak sector 
with gauge bosons, gaugino and higgsino, and the black points take into account the total 
corrections which include strong, weak and the contributions from the Higgs sector. As one 
can deduce from the plots, the corrections (top-stop, bottom-sbottom) coming from the 
strong interactions are independent of the triplet and singlet Higgs couplings, as expected, 
with a maximum split of 50 GeV respect to the tree-level mass eigenvalue. 

In the triplet-singlet extension we have four CP-even, three CP-odd neutral Higgs bosons 
and three charged Higgs bosons as shown in Eq. (3.7). These enhance both the Higgs and 
higgsino contributions to the radiative correction. The weak corrections (blue points) are 
dominated by the large number of higgsinos which contribute negatively to the mass and 
tend to increase for large values of the Higgs couplings (A t,s and K )- 

Finally, the black points show the sum of all the sectors, which are positive in sign, due 
to the large number of scalars contributing in the loop, with an extra factor of two for 
the charged Higgs bosons. This factor of two originates from the CW expression of the 
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Figure 6. The radiative corrections at one-loop for m, tl vs (a) A t, (b) A 5 and (c) k. The red points 
are only with strong (top-stop, bottom-sbottom) corrections, blue points are with weak corrections 
without the Higgs bosons (higgsinos, gauge boson and gauginos), (c) black points are the total 
(strong Tweak + Higgs bosons) corrections. 


potential, and accounts for their multiplicity (±). Such scalar contributions increase with 
the values of the corresponding couplings At, Ag, k. From Figure 6 one can immeditaley 
notice that the electroweak radiative corrections could be sufficient in order to fulfill the 
requirement of a ~ 125 GeV Higgs mass, without any contribution from the strong sector. 

To illustrate this point, in Figure 7(a) we have plotted the lightest CP-even neutral 
Higgs mass at one-loop versus the lighter stop mass (my,). We have used the same color 
coding conventions of the tree-level analysis. The red points are mostly doublets (> 90%), 
the green points are mostly triplet/singlet (> 90%) and blue points are mixed ones, as 
explained in section 3. The yellow band shows the Higgs mass range 123 < m/ ll < 127 
GeV. We notice that a ~ 125 GeV CP-even neutral Higgs could be achieved by requiring 
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Figure 7. The variation of the one-loop lightest CP-even Higgs mass with (a) the lightest 
stop mass rrif , and (b) with tan ,5, respectively. The yellow band shows the candidate Higgs mass 
123 < m hl < 127 GeV. 



a stop of very low mass, as low as 100 GeV. This is due to the presence of additional tree- 
level and radiative corrections from the Higgs sectors. Thus, in the case of extended SUSY 
scenarios like the TNMSSM, the discovery of a ~ 125 GeV Higgs boson does not put a 
stringent lower bound on the required SUSY mass scale, and one needs to rely on direct 
SUSY searches for that. 

In Figure 7(b) we present the dependency of the one-loop corrected Higgs mass of the 
lightest CP-even neutral Higgs on tan/3. The distribution of points is clearly concentrated 
at low values of tan/3 < 4. This is due to the additional contributions on the tree-level 
Higgs masses, which are maximal in the same region of tan/3 (see Eq. (3.8)). It is then clear 
that an extended Higgs sector reduces the amount of fine-tuning [12] needed in order to 
reproduce the mass of the discovered Higgs boson, compared to constrained supersymmetric 
scenarios. The latter, in general, require much larger supersymmetric mass scales beyond 
the few TeV [8] region. Compared to the pMSSM, this also represents an improvement, as 
it does not require large mixings in the stop masses in order to have the lighter stop mass 
below a TeV [9]. 

5.1 Hidden Higgs bosons 

Next we investigate the case in which we have one or more hidden Higgs bosons, lighter 
in mass than 125 GeV, scalars and/or pseudoscalars. In Figure 8 we present the mass 
correlations at one-loop for (a) m / ll — m^ 2 and (b) m ai — m a2 , where we have a CP-even 
candidate Higgs boson in the mass range 123 < m/,. < 127 GeV. The red points are mostly 
doublets (> 90%), green points are mostly triplets/singlets (> 90%) and blue points are 
mixed ones, as already explained. The green points have a high chance of evading the LEP 
bounds [6], showing that the possibility of having a hidden scalar sector is realistic, even 
after taking into account the radiative corrections to the mass spectrum. A closer inspection 
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Figure 8. The mass correlations at one-loop (a) — rrih 2 and (b) m ai — m a2 where we have a 

CP-even candidate Higgs in the mass range 123 < < 127 GeV. Red, blue and green points are 

defined as in Figure 2. 

of Figure 8 (a) reveals that there are points where both mt n and nih 2 are less than 100 GeV, 
showing that there is the possibility of having two CP-even hidden Higgs bosons. In that 
case /13 is the candidate Higgs of ~ 125 GeV. Similarly, Figure 8 (b) shows the possibility 
of having two hidden pseduoscalars. The arguments mentioned in section 3 will apply to 
the Higgs masses at one-loop as well. These imply that for < 123 GeV, the green 

points could evade the bounds from LEP and LHC, the red points would be ruled out and 
the blue points need to be carefully confronted with the data. In section 9 we analyse such 
scenarios in detail. The lightest pseudoscalar present in the spectrum, as we are going to 
discuss below, can play a significant role in cosmology. In fact, it is crucial in enhancing 
the dark matter annihilation cross-section, which is needed in order to get the correct dark 
matter relic in the universe [25]. 

6 /3-fuctions and the running of the couplings 

We have implemented the model in SARAH (version 4.5.5) [26] in order to generate the ver¬ 
tices and the model files for CalcHep [27], and generated the 0 functions for the dimension¬ 
less couplings and the other soft parameters. The 0 functions for A t,s,TS, n,9Y,9L,9c,Ut,b 
are given in the appendix A. 

To analyse the perturbativity of the couplings we have selected four different scenarios 
and identified the cut-off scale (A) in the renormalization group evolution, where one of the 
coupling hits the Landau pole and becomes non-perturbative (A*(A) = An). Figure 9(a) 
presents a mostly-triplet scenario at the electroweak scale as we choose A t = 0.8, A s,TS = 
0.1, k = 0.3 at the scale /to = M%. for tan 0 = 1.5 (solid lines) and tan 0 = 10 (dashed lines). 
For lower values of tan/3 (tan/3 = 1.5) the triplet coupling A t becomes non-perturbative 
already at scale of A ~ io 9 ~ 10 GeV, similarly to the behaviour shown in the triplet-extended 


16 - 








Ar = 0.8 


A 7 = A 5 = 0.8 




Ar = As = A ts = 0.8 


Ar = As = A ts = 0.8 k = 2.4 




In /i//i 0 


Figure 9. The running of the dimensionless Higgs couplings A t,s,ts an d k with the log of the ratio 
the scales (In/x/Ho) for tan/3 = 1.5 (solid lines) and tan/3 = 10 (dashed lines), where no = Mz- 
We have checked the corresponding variations for (a)A t = 0 .8, (b)Ar,s = 0.8, (c)A t,s,ts = 0.8 and 
(d)AT,s, = 0.8, k = 2.4 chosen at scale no respectively. 


MSSM [19, 28]. For larger values of tan /3 (tan /3 = 10) all the couplings remain perturbative 
up to the (Grand Unification) GUT scale (A ~ 10 16 GeV). 

Figure 9(b) presents the case where A t,s = 0.8 at no = Mz- We see that although 
the tan/3 dependency becomes less pronounced, the theory becomes non-perturbative at a 
relatively lower scale A 10 8 GeV. 

From Figure 9(c) it is evident that on top of A t and A s if we also choose A ts = 0.8 at 
Ho = Mz, the tan/3 dependency almost disappears. In this case the theory becomes more 
constrained with a cut-off scale A /'kj 10 6 GeV. 

Finally, Figure 9(d) illustrates the effect of a larger k value, the singlet self-coupling, 
with k = 2.4 at fiQ = Mz- The perturbative behaviour of the theory comes under question 
at a scale as low as 10 4 GeV. Such a large value of n at the electroweak scale thus restricts 
the upper scale of the theory to lay below 10 TeV, unless one extends the theory with an 
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extra sector 1 . Choosing relatively lower values of A ts and K would allow the theory to stay 
perturbative until 10 8 ~ 10 GeV even with A t,s as large as 0.8. The choice of larger values 
of A t,s increases the tree-level contributions to the Higgs mass (see Eq. (3.8)) as well as 
the radiative corrections, via the additional Higgs bosons exchanged in the loops. Both 
of these contributions reduce the amount of supersymmetric fine-tuning, assuming a Higgs 
boson of ~ 125 GeV in the spectrum, by a large amount, respect both to a normal and 
to a constrained MSSM scenario. Obviously, the addition of the triplet spoils the gauge 
coupling unification under the renormalization group evolution. This features is already 
evident in the triplet-extended MSSM [19, 28]. 


7 Fine-tuning 


The minimisation conditions given in Eq. (3.3) relate the Z boson mass to the soft breaking 
parameters in the form 


7 » /r2 2 2 

M z — /Uoft ~ /Vff 

1 


rn 2 H — tan 2 /3 m? H 


M = v s \ s - -^v T \ T , Aft = 2 tan2 p _ x 


ft is also convenient to introduce the additional parameter 

T = 


In ^sTft 

Msoft 


(7.1) 

(7.2) 


(7.3) 


characterizing the ratio between M| and /x 2 oft , which can be considered a measure of the 
fine-tuning. Unlike the MSSM, here the p e g parameter is generated spontaneously by 
the singlet and triplet vevs. Notice that while the triplet contribution is bounded by the p 
parameter [17], the singlet vev is unbounded and it may drive p, e s to a large value. Similarly, 
the soft parameters mn u ,H d , which are determined by the minimisation condition (3.3), can 
be very large, and thus can make /r 2 oft also large. Finally, to reproduce the Z boson mass 
we need large cancellations between these terms, which leads to the well know fine-tuning 
problem of the MSSM and of other supersymmetric scenarios. 

We show in Figure 10(a) plots of /U 2 oft and — /x 2 ff versus the singlet vev vs for tree-level 
candidate Higgs masses in the interval 120 < m ^ < 130 GeV. Figure 10(b) presents the 
same plots, but with m/ ll , the candidate Higgs mass, calculated at one-loop. The violet 
points represent /i 2 oft values for which A s,t > 0-5, and the points in blue refer to values 
of p 2 oit with A st < 0.5. The green points mark values of /r 2 ff with A $,t ^ 0.5, and the 
orange points refer to p? eS values with A s,T < 0.5. We see that for low A t,s both p 2 oft and 
—/r 2 ff (blue and orange points) are small, so that the required cancellation needed in order 
to reproduce the Z boson mass is also small. This leads to less fine-tuning, measured by 
T < 1. Unfortunately, such points are small in numbers in the tree-level case, since they 
require the extra contributions from the triplet and the singlet in order to reproduce the 

1 For the scan in Eq. 5.4 we select |k| < 3. The theoretical perturbativity of the parameter points have 

to be checked explicitly. 
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( a ) ( b ) 

Figure 10. The (a)tree-level and (b) one-loop level fine-tuning measures /x so ft and — /Xg ff versus the 
singlet vev vs for a candidate Higgs of mass between 120 < < 130 GeV respectively. The violet 

points represent values A s,t > 0.5 and blue points represent /Ltg oft values A s,t < 0.5. The green 

points represent values A s,t > 0.5 and the orange points represent values A s,t < 0.5. The 

red line shows the Z boson mass Mz- 


~ 125 GeV Higgs mass. For A t,s > 0.5 both and — (the violet and green points) are 
both very large, leading to large cancellations and thus to a fine-tuning parameter T ~ 5 
for “heff ~ 10 6 . 

Comparing Figure 10(a) and Figure 10(b) we see that the tree-level Higgs mass needs more 
fine-tuning as /Xg oft eff ~ 10 6 for large A t,s- The situation improves significantly at one-loop 
due to the contributions from the radiative corrections. This is due to the fact that there 
are more solutions with low values of A t, A t,s < 0-5, compared to tree-level and, on top of 
this, (for high and low A t,s) the required fine-tuning is reduced (J 7 < 2). This fine-tuning 
measure is a theoretical estimation but it is constrained from the lightest chargino mass 
bound from LEP (m~± > 104 GeV), which results in /r e fj > 104 GeV and T > 0.2. 

We have performed a run of Hi using the corresponding /3-functions for large Ar,s, 
from electroweak scale ( Mz ) up to a high-energy scale ~ 10 9,1 ° GeV, where the couplings 
become non-perturbative. It can be shown that Hd , /4ft no ^ blow up unless the 
couplings A t,s hit the Landau pole. The requirement of perturbativity of the evolution 
gives stronger bounds on the range of validity of the theory and the fine-tuning parameter 
is a good indicator at the electroweak scale. 

In the case of MSSM, the large effective quartic coupling comes from the storng SUSY 
sectors which also increase nrv\j and other parameters. However the situation changes in 
the case of extended Higgs sectors, which gives additional tree-level as well as quantum 
corrections to the Higgs masses. These reduce the demand for larger rrijj . In our case 
there is a singlet and a triplet which contribute largely at the tree-level for low tan /3 and 
also contribute at the quantum level. In the case of tree-level Higgs mass, the extra tree- 
level contributions demand very large A t,s ~ 0-8, which in turn make /i e fj very large and 
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so the fine-tuning. However in the case of Higgs mass at one-loop, the extra contributions 
from the extended Higgs sectors are shared by both tree-level and quantum corrections, 
which reduces the requirement of large At, s- This reduces /J> e g and so the fine-tuning T. 

8 A light pseudoscalar in the spectrum 

In the limit when the A parameters in Eq. (2.5) go to zero, the discrete Z 3 symmetry 
of the Lagrangian is promoted to a continuos 17(1) symmetry given by Eq. (8.1). This 
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Figure 11. The lightest CP-even Higgs boson mass mh x vs the lightest pseudo-scalar mass m ai at 
one-loop (top-stop and bottom-bottom corrections). The violet-yellow band presents the candidate 
Higgs mass 123 < < 127 GeV. The violet band specify the points with m ai < 1 MeV, where the 

a± —> e + e~ decay is kinematically forbidden. Red, blue and green points are defined as in Figure 2. 


symmetry is spontaneously broken by the vevs of the doublets, triplet and the singlet fields 
and should contain a physical massless pseudoscalar, a\, the Nambu-Goldstone boson of 
the symmetry. The soft breaking parameters will then lift the mass of ai, turning it into a 
pseudo-Goldstone mode whose mass will depend on the A j. The symmetry takes the form 

(H u , H d , f, S) ->• j*(H u , H d , f , S). (8.1) 

If this symmetry is softly broken by very small parameters A,- of 0(1 ) GeV, we get a very 
light pseudoscalar [10, 20] which could be investigated at cosmological level. Notice that 
the vector-like nature of the symmetry decouples this pseudoscalar from any anomalous 
behaviour. We are going to briefly investigate some features of the ai state in the context 
of the recent Higgs discovery and we will consider two different realizations. In the first 
case we consider a scenario where such continuous symmetry is broken very softly. In 
this case we choose the Ai parameters to be 0(1) GeV. We expect the pseudo-Goldstone 
boson to be very light, with a mass 0(1) GeV. In Figure 11 we show the mass correlation 
between the lightest CP-even neutral Higgs boson hi and the lightest massive CP-odd 
neutral Higgs boson a\. The red points are of doublet type, the green points represent 
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massive states of triplet/singlet type and the blue points represent the mixed contributions 
to the ai pseudoscalar. The violet-yellow band presents the region of parameter space 
where h\ is the candidate Higgs, with a mass 123 < m < 127 GeV. It is rather clear 
from Figure 11(a) that there plenty of points in the parameter space where there could be 
a hidden pseudoscalar Higgs boson along with or without a CP-even hidden scalar. Such a 
light pseudoscalar boson gets strong experimental bounds from LEP searches [6] and from 
the bottomonium decay rates [29]. Such light pseudoscalar in the mass range of 5.5 -14 
GeV, when it couples to fermions, gets strong bound from the recent CMS data at the LHC 
[30]. For triplet/singlet green points these bounds can be evaded quite easily since these 
states do not couple to gauge bosons (the Z boson in the case of a triplet) and to fermions. 
Of course, for real mass eigenstates the mixing between the doublet-triplet/singlet would 
be very crucial in the characterization of their allowed parameter space. 

For a mass of the ai of 0(100) MeV, the decay to 7r77, 7r7rvr could be an interesting 
channel to investigate in order to search for this state [31]. The simpler 2-particle decay 
channel ai —> nir is not allowed due to the CP conservation of the model. Due to the 
singlet/triplet mixing nature of this state, it decays into fermion pairs e + e _ , gp, rf, if 
kinematically allowed. Notice that there is no discrete symmetry to protect this state from 
decaying, preventing it from being a dark matter candidate [32], Now, if we choose the Ai 
parameters to be of 0(1) MeV then we get a very light pseudoscalar boson with mass of the 
same order, as shown in Figure 11(b). Such a bosons cannot decay to up, rr kinematically. 
Following the same reasoning, if its mass is < 1 MeV, then even the a\ —> e + W channel 
is not allowed and only the photon channel remains open to its decay. In this case the ai 
resembles an axion-like particle, and can be a dark matter candidate only if its lifetime 
is larger than the age of the universe [31, 33]. The pseudoscalar, in this case, couples to 
photons at one-loop, due to its doublet component which causes the state to have a direct 
interaction with the fermions. 

We recall that the effective lifetime of a light pseduoscalar decaying into two photons is 
given by Eq. (8.2) [33] 


647T 

9 2 ail rnl 


( 8 . 2 ) 


where g a , 77 is the effective pseduoscalar /fermion coupling which is proportional to the 
doublet-triplet/singlet mixing. Notice that the ai shares some of the behaviour of axion- 
like particles, which carry a mass that is unrelated to their typical decay constant, and as 
such are not described by a standard Peccei-Quinn construction. They find a consistent 
description in the context of extensions of the SM with extra anomalous abelian symmetries 
[34] [35] and carry a direct anomalous (contact) interaction to photons. Such interaction is 
absent in the case of a ai state. 

Along with the lightest neutralino of the TNMSSM, this particle can be a dark matter can¬ 
didate. In the supersymmetic context a similar scenario, with two dark matter candidates 
has been discussed in [36]. The role of this pseudoscalar state, in the context of the recent 
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results by FERMI about the 1-3 GeV excess gamma-ray signal from the galactic center [25] 
is under investigation for this model [37]. 


9 ~ 125 GeV Higgs and LHC data 

In this section we consider the one-loop Higgs mass spectrum, including only the correction 
coming from quarks and squarks, in light of recent results from the LHC [1-3] and the 
existing data from LEP [ 6 ]. In particular, we consider the uncertainties in the decay modes 
of the Higgs to WW*, ZZ* and 77 in a conservative way [1, 3]. We explore the scenario 
where one of the CP-even neutral scalars is the candidate ~ 125 GeV Higgs boson within 
the mass range 123 < m.y li < 127 GeV and investigate the possibilities of having one or 
more light scalars, CP-even and/or CP-odd, allowed by the LEP data and consistent with 
the recent Higgs decay branching fractions at the LHC. 

We just mention that in the TNMSSM the triplet and the singlet type Higgs bosons 
do not couple to the Z boson but the triplet couples to the W± bosons, which result in a 
modified hi W T vertices given by 

hi w ± = l - g\ {v u ni + v d ng + 4 v T nf A ), ( 9 . 1 ) 

where the rotation matrix Rfj is defined in Eq. (B.l). The vertices hi Z Z are given by 

hiZ Z = l -(g L cos 9 W + 9y sin 9 W f {v u fcfi + v d nf 2 ) , (9.2) 

where 9\y is the Weinberg angle. The Yukawa part of the superpotential is just the MSSM 
one. Hence the couplings of the CP-even sector to the up/down-type quarks and to the 
charged leptons are 


7 2 __ C* 

hiUU= -~^ VuRil ' 

(9.3) 

— Z c 

hi dd = --J=y d n i2 , 

(9.4) 

hi£t= 

(9.5) 


respectively. 

On the other hand, in the Higgs bosons decay into di-photons, there are more virtual 
particles which contribute in the loop compared to the SM. This is due to the enlarged 
Higgs and strong sectors which have a non-zero coupling with the photon. In particular 
there are three charginos tl iree charged Higgs bosons (Kf 23 ), the stops (£ 72 ) and 

the sbottoms ( 61 , 2 ). Compared to the MSSM and the NMSSM we have two additional 
charged Higgs bosons and one additional chargino which contribute to the decay. The 
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decay rate in the di-photon channel is given by 


F (h ->• 77 ) = 
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where JVj s are the color number of fermion and scalars, Qf t s are the electric charges, in 

TTlP 1 

unit of |e|, of the fermions and scalars, and A)> ^ 1/2 an d A\ are the spin-0, 

spin-l /2 and spin -1 loop functions 

4)0*0 = — 4 ( x ~ /( x )) > ( 9 - 7 ) 

x z 

a i/ 2 {x) = ^(x +(x-l)f(x)) , (9.8) 

4(x) = — (2 x 2 + 3 x + 3(2 x — 1 )/(x)) , (9.9) 

with the analytic continuations 
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In the limit of heavy particles in the loop, we have Aq —> 1/3, A i / 2 —> 4/3 and A] —>• —7. 
Using the expression above, we study the discovered Higgs boson (h. 125 ) decay rate to 
di-photon in this model. We also check the consistency of light scalar(s) and/or light 
pseudoscalar(s) with the current data at the LHC and the older LEP data. Such analysis is 
presented in Figure 12. Figure 12(a) shows such hidden Higgs scenarios with one a\ and/or 
one h\ below 123 GeV, which find significant realizations. We first consider the results 
coming from both CMS and ATLAS in the decay of the Higgs to WW* and ZZ* modes at 
lcr [1-3] and also consider the cross-section bounds from LEP [ 6 ]. The allowed mass values 
are shown as red points for which the lightest CP-even Higgs boson (hi) is the detected 
Higgs at ~ 125 GeV. Cleary we see that there are many light pseudo-scalars (< 100 GeV) 
which are allowed. The orange points present the scenario where mh 2 ~ mi 25 and which 
leaves both h\ and ai hidden (< 125 GeV). 

We have performed additional tests of such points and compared them with the results 
from the decay of the Higgs boson to di-photons at the LHC, both from CMS [38] and 
ATLAS [39]. The red points (with one hidden Higgs boson) which satisfy /1125 —> 77 at 
lcr level, are marked as green points. The orange points (with two hidden Higgs bosons) 
when allowed at 1 <j level, have been marked as blue points. Notice that all the points in 
Figure 12 are allowed at lcr by the WW *, ZZ* channels and at 2(7 by the 77 channel. These 
requirements automatically brings the fermionic decay modes closer to the SM expectation. 
Of course the uncertainties of these decay widths give us a room for h\ 25 — > a\ai/h\h\. 

Notice also the presence of a very light pseudoscalar mass values near a\ ~ 0. Fig¬ 
ure 12(b) is a zoom of this region, where such solutions are shown for m ai < 1 MeV. The 
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Figure 12. The lightest CP-even Higgs boson mass vs the lightest pseudo-scalar mass m ai 
at one-loop (top-stop and bottom-bottom corrections) consistent with the Higgs data from CMS, 
ATLAS and LEP. The red points corresponds to the case where rrp (1 ~ TO 125 , the orange points 
correspond to mass values of and m ai where rrih 2 ~ 771125 and all of them satisfy the ZZ *, 
WW* bounds at la and 77 bound at 2er level from both CMS and ATLAS. The red (orange) points 
which satisfy the 77 result at la are marked green (blue). Very light pseudoscalar masses m ai < 1 
MeV are shown in panel (b), which is a zoom of the small mass region of (a). 


points in this case correspond to possible ai states which do not decay into any charged 
fermion pair (m ai < 2 m e ) and have an interesting phenomenology, as briefly pointed out 
in section 8 . The fact that such mass values only allow a decay of this particle to two 
photons via doublet mixing mediated by a fermion loop, makes the a\ a possible dark mat¬ 
ter candidate, being long lived. Two hidden Higgs bosons render the phenomenology very 
interesting, allowing both the h\ 2 § —> aiai and the /1125 — > h\h\ decay channels [40, 41]. In 
Figure 13 we show some of the points in this model as benchmark points (BMP’s), which 
are allowed both by LHC [1, 3] and LEP [ 6 ] data. The neutral Higgs spectrum has been 
calculated at one-loop order and the remaining states at tree-level. Figure 13(a) shows a 
point (BMP1) where we have a hidden pseudoscalar (ai) with mass 0(1W 1 ) MeV and an¬ 
other triplet/singlet-like hidden CP-even scalar (h\) with a mass around ~ 93 GeV. In this 
case the candidate Higgs boson is / 12 , taken around 125 GeV. This point also have a triplet 
type very light charged Higgs boson at a mass around 90 GeV, which is not excluded by the 
recent charged Higgs bounds from the LHC [22], Figure 13(b) shows a benchmark point 
(BMP2) where we have a pseudoscalar around 37 GeV, and the lightest scalar and charged 
Higgs bosons around 100 GeV. Figure 13(c) shows a trivial (SM-like) solution where we 
have a doublet-type CP-even Higgs around ~ 125 GeV, with the other states decoupled. In 
the next study we are going to analyse such points through a detailed collider simulation 

[371 ■ 
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Figure 13. We show the benchmark points of the model which are allowed both by LHC [1, 3] 
and LEP[6] data. The neutral Higgs spectrum has been calculated at one-loop and the rest of the 
spectrum at tree-level. 


10 Phenomenology of the TNMSSM 

The TNMSSM extends the Higgs sector as well as the electroweak chargino-neutralino 
sectors by additional higgsino contributions. Both the triplet and singlet fields do not 
couple to the fermions but affect the phenomenology to a large extent. In the context of 
the recent Higgs discovery, searches for additional Higgs bosons, both neutral and charged, 
are timely. In particular, if an extended Higgs sector will be discovered at the LHC, it will be 
crucial to determine the gauge representation which such states belong to, by investigating 
its allowed decays modes. 

We have seen from Eq. (2.2), that af = 0 hypercharged triplet couples to the 
bosons and contributes to their mass. On the other hand, the singlet does not directly 
couple to any of the gauge bosons. In the case of Higgs mass eigenstates which carry a 
doublet-triplet-singlet mixing, we need to look either for their direct production processes 
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Figure 14. The new and modified production channels for the Higgs bosons at the LHC. 


at the LHC or take into consideration the possibility of their cascade production from other 
Higgses or supersymmetric particles. 

10.1 Productions 

We have detailed a model with a rich Higgs sector with additional Higgs bosons of triplet 
and singlet type. We recall that the relevant production processes of a Higgs boson which 
is a SU(2) doublet at the LHC [42, 43] are the gluon-gluon fusion (GGF) and vector boson 
fusion channels, followed by the channels of associated production of gauge bosons and 
fermions. In our case, the production channels for the new Higgs bosons are different, due 
to their different couplings to the gauge bosons and fermions. We list below the possible 
additional production channels for the neutral and charged Higgs bosons at the LHC. 

• Neutral Higgs boson production in association with charged Higgs boson: The triplet 
only couples to boson. Thus a neutral Higgs (doublet or triplet) can be produced 
in association with a charged Higgs boson (doublet or triplet) via a W± exchange. 
As shown in Figure 14(a) a light in mass and charged Higgs boson in the TNMSSM 
can be easily explored by this production channel qq' —> h t . 

• Neutral Higgs boson production in association with W^: A triplet or a doublet type 

neutral Higgs boson can be produced via qq' —>• W±hi as shown in Figure 14(b). A 
triplet admixture modifies the hi — — W T couplings by an additional term propor¬ 

tional to the vev of the triplet. 

• Charged Higgs boson production in association with W±\ Triplet of Y = 0, ±2 hy¬ 
percharge has a non-zero tree-level coupling to Z — W± — hf. This leads to additional 
contributions to qq —> W^tif as shown in Figure 14(c). 
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• Production of charged Higgs boson in vector boson fusion: The non-zero Z—W^ — hf 
coupling leads to vector boson fusion (Z,W fusion) which produces a charged Higgs 
boson as shown in Figure 14(d). This mode is absent in 2-Higgs doublet models 
(2HDM), in the MSSM and in the NMSSM. This is a unique feature of the Y = 0, ±2 
hypercharge, triplet-extended scenarios. 

• Singlet Higgs production: The singlet in this model is not charged under any of the 
gauge groups, and hence the direct production of such a singlet at the LHC is impos¬ 
sible. Gauging this additional singlet with the inclusion of an extra additional U{1)' 
gauge group would open new production channels via the additional gauge boson 
(Z 1 ). Most of the extra Z' models get a bound on the Z' mass, mz> > 2.79 TeV [44] 
which makes such channels less promising at the LHC. In our case such a singlet type 
Higgs boson would only be produced via mixing with the Higgs bosons of doublet and 
triplet type. 

10.2 Decays 

The smoking gun signatures for the model would be the decays of the doublet, triplet and 
singlet states that are produced. Different F-term contributions can generate these types 
of mixing and corresponding decay vertices. We list all the vertices in the gauge basis in 
appendix B. The vertices in the mass eigenstate basis can be found by the rotations given 
in Eq. (B.l). 



(c) (d) (e) 

Figure 15. The new and modified decay channels of the Higgs bosons at the LHC. 


• Higgs decays to Higgs pairs: The candidate Higgs around the 125 GeV mass in this 
case can decay into two hidden Higgs bosons if this channel is kinematically allowed 
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as can be seen in Figure 15(a). Such hidden Higgs boson(s) could be both scalar and 
pseudoscalar in nature. The discovered Higgs is however 99% CP-even [2], which for¬ 
bids any CP-violating decay of the nature h\2b —> aihj. However, the CP-conserving 
decays like h ±25 —> cuaj and/or h\ 2 5 —> h{hj are allowed. Such decays should be 
carefully investigated on the basis of the current Higgs data at the LHC. If the two 
light Higgs bosons are mostly singlet or triplet, then it is easy to evade the bounds 
from LEP [ 6 ]. As we have already pointed out, a singlet Higgs does not couple to 
any of SM gauge bosons and even the triplet type does not couple to the Z boson. 
Such a light Higgs boson could decay into r pairs only through the mixing with the 
doublet type Higgs bosons, since neither the singlet nor the triplet couple to fermions 
(see Eq. (2.2)). The mixing angle is also constrained by data on bottomonium decay, 
for a very light neutral Higgs boson (< 8 GeV) [29]. 

The decay of a Higgs boson into other Higgs bosons depends on the cubic coupling, 
which is proportional to the vevs of Higgs fields, and thus it is very sensitive to the 
values of Uj. It is therefore requires an analysis of the allowed decay widths of the 
Higgs boson into a light Higgs pair using LHC data [40]. 

• Higgs decays to W±W T '. The triplet couples to W± via its non-zero SU(2) charge 
which is at variance respect to the analogous coupling of the doublet, as can be seen 
in Eq. B.2. This will modify the decay width of hi —> WW (Figure 15(b)). The recent 
data show that there is some disagreement and uncertainties between the CMS [1, 2 ] 
and ATLAS [3] results in the /1125 — > WW* channels. The measurement of this decay 
channel thus becomes even more crucial under the assumption of a triplet mixture. 

• Charged Higgs decays to ZW±: We know that the triplet type charged Higgs has 
a non-zero tree-level coupling to Z W, for a non-zero triplet vev, as shown in Fig¬ 
ure 15(c). This opens up the possible decay modes hf —> ZW±, which are absent in 
the 2HDM and in the MSSM at tree-level. 

• Charged Higgs decays to hj(aj)W±: A doublet or triplet type Higgs boson can decay 
to a lighter neutral Higgs and a W± (Figure 15(d)). A possibility of a very light 
triplet-singlet-like neutral Higgs makes this decay mode more interesting compared 
to the case of the CP-violating MSSM [45]. 

• Higgs decays to ajZ: In the MSSM the odd and heavy Higgs bosons are almost de¬ 
generate, so hi —>• a.jZ is not kinematically allowed. The introduction of a triplet and 
of a singlet adds two more massive CP-odd Higgs bosons, and the degeneracy is lifted. 
In this case we have a relatively lighter CP-odd Higgs state a* which makes hi —> ajZ 
possible, as shown in Figure 15(e). This scenarios is also possible in the context of 
the CP-violating MSSM, where we have a very light pseudoscalar Higgs boson due 
to the large mixing between the Higgs CP eigenstates [46] and in the NMSSM, for 
having an additional scalar [ 10 ]. 

• Higgs decays to fermion pairs: In a scenario where a triplet or/and singlet type Higgs 
boson decays to gauge bosons and other Higgses are kinematically forbidden, the only 
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permitted decays are into light fermion pairs, viz, bb, rr and fjfi. Even such decays 
are only possible by a mixing with doublet type Higgs bosons. When such mixing 
angles are very small this can results into some displaced charged leptonic signatures. 

10.3 Possible signatures 

The unusual production and decay channels lead to some really interesting phenomenology 
which could be tested in the next run of the LHC and at future colliders. From the testability 
point of view, one could use the data form the discovered Higgs boson ~ 125 GeV in order 
to get bounds from the Higgs decaying to Higgs boson pair [40], or the existing bounds 
from LEP [6] for two Higgs bosons productions. We have already taken into account these 
bounds by ensuring that the hidden Higgs boson is mostly of singlet or of triplet type. 
Given the uncertainty in the Higgs decay branching fractions in different modes and the 
absence of direct bounds on the non-standard decays of Higgs boson to Higgs boson pair 
(h \25 —>■ aidj/hihj ), this remains phenomenologically an interesting scenario. Below we list 
different possible signatures that could be tested in the LHC with 13/14 TeV. 

• The singlet and doublet F-terms generate the doublet-triplet-triplet vertex which is 
proportional to XsXts and Ay. This would provide a signature of a doublet type 
Higgs decaying into two triplet type Higgs bosons, which, in turn, do not decay into 
fermions. Similarly the F-terms of H u and H 4 generate vertices involving triplet- 
singlet-doublet which are proportional to XtXs ■ The F-term of triplet type also 
contributes to this mixing, which is proportional to XtXts • Thus the relative sign 
between the two contributions become important. The vertex is given in Eq. B.6. 

In the case of a ~ 125 GeV Higgs boson, this can decay into two triplet-like hidden 
scalars or pseudoscalars, which in turn decay into off-shell VF^s only. This type of 
decays can be looked for by searching for very soft jets or leptons coming from the 
off-shell VF^s. The signatures could be the 4£+ fir or 4 j + 21+ fix channels, where 
the jets and the leptons are very soft. On the other hand, both the triplet and the 
singlet hidden Higgses can decay to fermion pairs (bb, cc, e + e~, fj.fi, rf) via the mixing 
with doublets. The recent bounds on these non-standard decays has been calculated 
for the LHC [41]. Such decays give 4 1, 2b+ 21 final states, where the leptons are very 
soft. For the triplet type hidden Higgs bosons it would be interesting to analyze the 
competition between the four-body and the two-body decays (which depend on the 
triplet-doublet mixing). Demanding for the presence of softer leptons and jets in the 
final states, allows to reduce the SM backgrounds at the LHC. If the mixing is very 
small, this could lead to displaced charged leptonic final states, similar to those of a 
Higgs boson decay in a li-parity violating supersymmetric scenario [47]. Due to the 
coupling both with the up and the down type doublets, this coupling could be tested 
both at a low and a high tan/3. 

• The singlet does not contribute to charged Higgs mass eigenstates, so the charged 
Higgs bosons could be either of doublet or triplet nature. In the case of a heavy 
doublet type, the heavier charged Higgs can decay to a triplet type a neutral Higgs 
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(CP even or odd) and a triplet type charged Higgs {H± d —> T 0 Tj ± 2 ) ( see Appendix B). 
The coupling is proportional to — A j,). The lighter triplet type charged Higgs then 
mostly decays into on-shell or off-shell Z W±. This is a generic signature for Y = 0, ±2 
hypercharge triplets with non-zero triplet vev, which breaks the custodial symmetry 
of the Higgs potential. The relatively lighter triplet (either CP-odd or even) neutral 
Higgs can decay via an on/off-shell W± boson pair, which leads to leptonic final 
states. The final states with multi-lepton(> 3£), multi-jet(> 4) and missing energy, 
could be the signature for this model. Depending on the off-shell decays, few leptons 
or a jet could be softer in energy. 

• In other cases a triplet type heavier charged Higgs can decay into a doublet type 
neutral Higgs and a triplet type charged Higgs. These couplings are proportional to 
(^p — Af ) and can give rise to 3 i + 26+ pr and 3£ + 2r+ pr final states. Here the 6 
and r pairs expected from the neutral doublet type Higgs boson decay. 

• Unlike to the neutral Higgs bosons, the up and down type charged Higgs bosons 
doublet only mix with the triplets. The couplings are again proportional to a combi¬ 
nation of Xs^t an d XtXts • In this case the doublet (triplet) charged Higgs state will 
decay into a triplet (doublet) charged Higgs and a singlet neutral Higgs boson. As 
the singlet is not coupled to any SM particles, it can only decay through mixing with 
doublets and triplets. Decays of such singlets to leptons (in the case of mixing with 
doublets) and off-shell or on-shell VF^-pair will be determined by the mixing only. In 
a fine-tuned region where such mixing is very low this decay channel can lead to a 
displaced vertex of charged leptons, whose measurement can give information about 
such a mixing. 

11 Discussions and Conclusions 

In this paper we have considered a scenario with an extended Higgs sector characterized 
by a Y = 0 hypercharge SU( 2) triplet and a gauge singlet superfields, along with the 
remaining MSSM superfields. The triplet vev is restricted by the p parameter, hence the 
g e ff is generated spontaneously mostly by the singlet vevs. In models with gauged U(l)' 
symmetry the singlet could be invoked in the mass generation of the extra gauge boson Z l 
by spontaneous symmetry breaking. This would require a large singlet vev vs, due to the 
recent bounds on extra Z' coming from the analysis at the LHC [44], 

We have first investigated the masses of the Higgs sector of the model at tree-level. 
The lightest tree-level Higgs state, in this case, is not bounded to lay below Mz, due to the 
additional contributions from the triplet and the singlet, which are proportional to their 
respective couplings and are enhanced at low tan/3. This allows to reduce the size of the 
quantum correction needed in order to reach the ~ 125 GeV at one-loop, compared to the 
MSSM or to others constrained MSSM scenarios. Then we have extended our analysis at 
one-loop level. The one-loop Higgs with mass around ~ 125 GeV puts some indirect bounds 
on the masses of the particles contributing in the radiative corrections. For this purpose we 
have included the one-loop contributions using the Coleman-Weinberg potential. We have 
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also presented results for the neutralino, and charginos spectra, together with the stop and 
sbottom mass matrices. We have calculated full one-loop Higgs masses considering all the 
weak sectors and the strong sectors. We also showed that the gauge boson-gaugino-higgsino 
sectors mostly contribute negatively to the mass eigenstates, while the stop-top, sbottom- 
bottom and Higgs sectors contribute positively. Due to the large number of scalars, seven 
neutral and three charged Higgs bosons, the Higgs self corrections can be larger than the 
strong corrections in the large A t,s limit. This substantially reduces the indirect lower 
bounds on the stop and sbottom masses. Thus in TNMSSM the discovery of a ~ 125 GeV 
Higgs boson does not put a stringent lower bound on the stop and sbottom masses, and 
one has to rely on direct search results for the lower bounds on the SUSY mass scale. 

We have implemented the model in SARAH3.5 [26] in order to generate the vertices and 
other model hies for CalcHEP [27]. The beta-functions have also been generated at one- 
loop. We have addressed the issues of perturbativity of the couplings at the higher scale, as 
we have run the corresponding renormalization group equations from the electroweak scale 
up. This has shown that the couplings of the model at the electroweak scale need to be 
restricted to certain values. For example, even with a value of A t,s ~ 0.8 at the electroweak 
scale, the theory remains perturbative up to 10 8 ~ 10 GeV. Setting all the couplings at a value 
(A ts ~ 0.8, k ~ 2.4) the upper scale in the perturbative evolution gets lowered to 10 4 ~ 6 
GeV. The issue of fine-tuning at the electroweak scale has been discussed in this context. 
We have seen that although the tree-level mass spectrum is highly fine-tuned for larger 
A t,Si the amount of fine tuning is reduced after the inclusion of the radiative corrections. 

The prospects for hidden Higgs(es), which are scalars and/or pseudoscalars of mass 
lower than the current Higgs mass, has been discussed quite thoroughly. We have seen 
that in the rich Higgs spectrum of the model there are several possibilities for having one 
or more hidden neutral Higgs bosons (< 125 GeV) both CP-even and CP-odd. A special 
scenario emerges when we break the continuous U( 1) symmetry softly by the parameters 
Ai. This leads to the appearance of a very light pseudoscalar state of 0(1 ) GeV to 0(1) 
MeV in mass, which has its own interesting phenomenology. 

Finally, we have discussed the doublet-triplet-singlet mixing which influences the pro¬ 
ductions and decays of neutral and charged Higgs bosons at the LHC. The existence of a 
hf — W^ — Z tree-level vertex, due to the triplet, impacts both the production as well as the 
decay channels of the charged Higgs bosons [15]. In the presence of a light pseudoscalar, the 
hi Zcij channel is a possibility due to the very light mass of the pseudoscalar (s). Both 
the triplet and the singlet states do not couple to the fermions, which leads to some very 
interesting phenomenology. This property also has an impact on rare decays like b —> /z/z 
and b —> s'y [12, 31]. Given the rich phenomenology and the specific predictions of this 
model, the current analysis at the LHC and future colliders could be able to test and shed 
a light on this scenario by looking at its interesting signatures. 

Acknowledgements 

P.B. wants to thank Eung Jin Chun and Ajit Mohan Shrivastava for the discussion regarding 
light pseudo-scalar as hot or cold dark matter candidate. The work of C.C. is supported 


- 31 


by a The Leverhulme Trust Visiting Professorship at the University of Southampton at the 
STAG Research Centre and Mathematical Sciences. He thanks Kostas Skenderis and all 
the member of the Centre for the kind hospitality. 


A RG equations 


We list the RG equations at one-loop for the dimensionless coupling A t,S,TS, 9y,9l- 

Here t = In . where fj, is the running scale and /jq is the initial scale. 
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B Higgs Vertices 


The Higgs boson vertices are given in the gauge basis and the vertices in the Higgs boson 
mass basis can be obtained by the rotation defined in in Eq. (B.l). 7 Zfj is the rotation matrix 
for CP-even neutral Higgs boson and H, L = (H® r , H® r , S r , TJ?), hi = (hi, / 12 , / 13 , h. 4 ) are the 
CP-even neutral Higgs bosons in the gauge and mass basis respectively. For the pseudo- 
scalar we used the rotation matrix 1Z P which rotates the CP-odd neutral Higgs boson from 
their gauge basis, Ai = (H® -, H^ i: S 4 , Tf ) to the mass basis a* = ( 00 , 01 , 02 , 03 )- T^ C is 
the corresponding rotation matrix that rotates the gauge basis H p = (_ff+, T^~, H7, T^~) 
to mass eigenstates hf = (hg , h±, h ±, hg ). Here ao and hg are the neutral and charged 
Goldstone bosons which give masses to the Z and bosons respectively. 


hi = KfjHj 

O* = 'R'ijAj 

hf = n%Hf 

(B.l) 

T°H° ,H° f At + XtsX t v s XsXtsvt • 

± r ±± u,r/i 1 ± d,r/i ■ 2 ~ ^ l - > ^ ’ 

(B.2) 

^0 zjO u-o . x tsX T v s A t x , .. . 

1 i H u,r/i H d,i/r ■ ^2 2 X TSX S VT, 

(B.3) 

zjO u 0 . \ _ f V T X T V S X 

**u/d, r/i **u/d,r/i ' * T V 4 2\/2 ’ 

(B.4) 

0^0 7_t0 ttO . n. 

n u/d,r n u/d,i • U ’ 

(B.5) 

H u/d,r ■ -p^(v d / u XTS ~ V u / d Xs)XT J 

(B.6) 

T r H u/d,i S i : p^ v u/dXTXTS] 

(B.7) 

-^i H u /d,r/i Si/r ■ ~~ i V u/dXs T v d/uXTs)XT , 

(B.8) 

T® Sr Sr : VT Xts(k + 2 Ats); 

(B.9) 

T° S r Si : 2 k,Xtsvt‘, 

(B.10) 

T r ° Sj Si : nr Ats( 2 Ats - «); 

(B.ll) 

T r T r S r : A T s u s (k + 2A r s) + ; 

(B.12) 

T r ° Tj° Si : 2 k us Ats - V 2 A T s\ 

(B.13) 

Ti T t Sr : Ats ^s ( 2 Ats «) ^ ’ 

(B.14) 

7-1O 7-1O ry-iO \2 

1 r 1 r/i 1 r/i ■ V T *TSi 

(B.15) 

S r S r S r : A * + k 2 v s - 
3^2 

(B.16) 


33 - 









S r Si Si : k 2 v s -^|; 

(B.17) 

c uO uO . \2 VtXt^S. 

°r n u/d,r/i n u/d,r/i ' 2 ^ SA 5 5 

(B.18) 

e yo Tj-o . v t XtXs . -4s. 

**u/d,r/i **d/u,r/i ’ y^2’ 

(B.19) 

c rrO As . , V T X T X T S 

^i n u/d,r/i n d/u,i/r ■ ^ K,AsVs + ^ ’ 

(B.20) 

SiH° u/d , r H° /dji : 0 ; 

(B.21) 


H u/d,r H u/d,r/i H u/d,r/i : g (Sy + 9w) v u/d 

(B.22) 

H u/d,r H d/u,r/i H d/u,r/i : 2 (A! + 9 -^T “ ^(Sy + 9w)) v u/d', 

(B.23) 

uO ttO ttO . rv. 

**u/d,i ^u/d,r **d/u,i ' 

(B.24) 

H u/d,r T r T r : -^ u / d >^ - ~V d / u X S XTS', 

(B.25) 

H u/d,r T i T i '■ ^ v u/dXr + 2 V d/ U XsX T s', 

(B.26) 


H u/d,i T r T i : -Vd/uXsXrs; 

(B.27) 

Hu/d,r A- ’ 2 ^s( v u/dXs ~ v d/u K )'i 

(B.28) 

^u/d,r A A • 2 Xsi^u/dXs 4“ V d j u K), 

(B.29) 

H u/d,i S rSi : -v d/u nX s \ 

(B.30) 

- 34- 




T«H+HJ 

^ 1 > 

(B.31) 


: A ts(vs\t + v^vtAs); 

(B.32) 

^<7 (<40 

A t , A tvt , 

: 2 {XsVS + ^2 ); 

(B.33) 


: 0; 

(B.34) 

r»r-/ + (!■-/+)* 

Sl 

: 7T t; 

(B.35) 

T° Tf T+ 

: 0; 

(B.36) 

T° (Tf)* (Tf )' 

: V2vt{2\ts ~ 9l)\ 

(B.37) 

r %2 + cA/~)' 

: 0; 

(B.38) 

rpO rp — / + rp— / + 

1 J 1/2 M/2 

: 0; 

(B.39) 

rpO rp~ /+ pf + / — 

1 1 l/2 n u/d 

: 0; 

(B.40) 

T° (Ti"/ + )t «/„")* 

■ T Vu ' d M gly, 

(B.41) 

r %2 + (</") t 

: 0; 

(B.42) 

s KICK'D' 

: a s (a t ^ + 

(B.43) 

sif/r (K/;)' 

: -p.\s^ T v u / d ; 

(B.44) 


References 

[1] G. Aad et al. [ATLAS and CMS Collaborations], arXiv: 1503.07589 [lrep-ex]. S. Chatrchyan et 
al. [CMS Collaboration], JHEP 1401 (2014) 096 [arXiv:1312.1129 [hep-ex]]. S. Chatrchyan et 
al. [CMS Collaboration], Phys. Rev. D 89 (2014) 9, 092007 [arXiv: 1312.5353 [lrep-ex]]. 

[2] V. Khachatryan et al. [CMS Collaboration], arXiv: 1412.8662 [lrep-ex], [CMS Collaboration], 
CMS-PAS-HIG-13-002. 

[3] ATLAS-CONF-2015-007 G. Aad et al. [ATLAS Collaboration], Phys. Rev. D 91 (2015) 1, 
012006 [arXiv: 1408.5191 [lrep-ex]]. G. Aad et al. (ATLAS Collaboration], arXiv: 1412.2641 
[lrep-ex]. G. Aad et al. [ATLAS Collaboration], Phys. Rev. D 90 (2014) 5, 052004 
[arXiv:1406.3827 [hep-ex]]. 

[4] E. Witten, Nucl. Phys. B 188 (1981) 513. S. Dimopoulos and H. Georgi, Nucl. Phys. B 193 
(1981) 150. E. Witten, Phys. Lett. B 105 (1981) 267. R. K. Kaul and P. Majumdar, Nucl. 
Phys. B 199 (1982) 36. N. Sakai, Z. Phys. C 11 (1981) 153. 

[5] A. H. Chamseddine, R. L. Arnowitt and P. Nath, Phys. Rev. Lett. 49 (1982) 970. G. L. Kane, 
C. F. Kolda, L. Roszkowski and J. D. Wells, Phys. Rev. D 49 (1994) 6173 [hep-ph/9312272], 

[6] R. Barate et al. [LEP Working Group for Higgs boson searches and ALEPH and DELPHI 
and L3 and OPAL Collaborations], Phys. Lett. B 565 (2003) 61 [hep-ex/0306033]. S. Schael 
et al. [ALEPH and DELPHI and L3 and OPAL and LEP Working Group for Higgs Boson 
Searches Collaborations], Eur. Phys. J. C 47 (2006) 547 [lrep-ex/0602042]. 


35 - 



[7] R. Barbier, C. Berat, M. Besancon, M. Chemtob, A. Deandrea, E. Dudas, P. Fayet and 

S. Lavignac, G. Moreau, E. Perez, Y. Sirois Phys. Rept. 420 (2005) 1 [hep-ph/0406039]. 

[8] H. Baer, V. Barger and A. Mustafayev, JHEP 1205 (2012) 091 [arXiv: 1202.4038 [lrep-ph]]. 

J. Ellis and K. A. Olive, Eur. Phys. J. C 72 (2012) 2005 [arXiv:1202.3262 [lrep-ph]]. P. Nath, 
Int. J. Mod. Phys. A 27 (2012) 1230029 [arXiv:1210.0520 [hep-ph]]. 

[9] M. Carena, S. Gori, N. R. Shah and C. E. M. Wagner, JHEP 1203 (2012) 014 
[arXiv:1112.3336 [hep-ph]]. L. J. Hall, D. Pinner and J. T. Ruderman, JHEP 1204 (2012) 

131 [arXiv:1112.2703 [hep-ph]]. A. Arbey, M. Battaglia, A. Djouadi, F. Mahmoudi and 

J. Quevillon, Phys. Lett. B 708 (2012) 162 [arXiv:1112.3028 [hep-ph]]. 

[10] N. E. Bomark, S. Moretti, S. Munir and L. Roszkowski, JHEP 1502 (2015) 044 
[arXiv:1409.8393 [hep-ph]]. G. G. Ross and K. Schmidt-Hoberg, Nucl. Phys. B 862 (2012) 
710 [arXiv:1108.1284 [hep-ph]]. K. S. Jeong, Y. Shoji and M. Yamaguchi, JHEP 1204 (2012) 
022 [arXiv:1112.1014 [lrep-ph]]. U. Ellwanger, Eur. Phys. J. C 71 (2011) 1782 
[arXiv:1108.0157 [hep-ph]]. U. Ellwanger, C. Hugonie and A. M. Teixeira, Phys. Rept. 496 
(2010) 1 [arXiv:0910.1785 [hep-ph]]. 

[11] J. R. Espinosa and M. Quiros, Nucl. Phys. B 384 (1992) 113. J. R. Espinosa and M. Quiros, 
Phys. Lett. B 279 (1992) 92. 

A. Delgado, M. Garcia-Pepin, B. Ostdiek and M. Quiros, arXiv: 1504.02486 [hep-ph], L. Cort, 
M. Garcia and M. Quiros, Phys. Rev. D 88 (2013) 7, 075010 [arXiv: 1308.4025 [hep-ph]]. 

M. Frank, D. K. Ghosh, K. Huitu, S. K. Rai, I. Saha and H. Waltari, Phys. Rev. D 90 (2014) 
11, 115021 [arXiv: 1408.2423 [lrep-ph]]. M. Garcia-Pepin, S. Gori, M. Quiros, R. Vega, 

R. Vega-Morales and T. T. Yu, Phys. Rev. D 91 (2015) 1, 015016 [arXiv: 1409.5737 [lrep-ph]]. 
M. Das, S. Di Chiara and S. Roy, Phys. Rev. D 91 (2015) 5, 055013 [arXiv:1501.02667 
[hep-ph]]. 

[12] P. Bandyopadhyay, S. Di Chiara, K. Huitu and A. S. Kei£jeli, JHEP 1411 (2014) 062 
[arXiv:1407.4836 [hep-ph]]. P. Bandyopadhyay, K. Huitu and A. Sabanci, JHEP 1310 (2013) 
091 [arXiv: 1306.4530 [hep-ph]]. 

[13] S. Di Chiara and K. Hsieh, Phys. Rev. D 78 (2008) 055016 [arXiv:0805.2623 [hep-ph]]. 

[14] J. L. Diaz-Cruz, J. Hernandez-Sanchez, S. Moretti and A. Rosado, Phys. Rev. D 77 (2008) 
035007 [arXiv:0710.4169 [lrep-ph]]. 

[15] P. Bandyopadhyay, K. Huitu and A. S. Keceli, JHEP 1505 (2015) 026 [arXiv:1412.7359 
[hep-ph]]. E. Barradas-Guevara, O. Felix-Beltran and A. Rosado, Phys. Rev. D 71 (2005) 
073004 [hep-ph/0408196]. 

[16] K. rn. Cheung, R. J. N. Phillips and A. Pilaftsis, Phys. Rev. D 51 (1995) 4731 
[hep-ph/9411333]. R. Godbole, B. Mukhopadhyaya and M. Nowakowski, Phys. Lett. B 352 
(1995) 388 [hep-ph/9411324]. K. Cheung and D. K. Ghosh, JHEP 0211 (2002) 048 
[hep-ph /0208254]. 

[17] J. Beringer et al. [Particle Data Group Collaboration], Phys. Rev. D 86 (2012) 010001. 

[18] T. Basak and S. Mohanty, Phys. Rev. D 86, 075031 (2012) [arXiv:1204.6592 [hep-ph]]. 

T. Basak and S. Mohanty, JHEP 1308 (2013) 020 [arXiv:1304.6856 [hep-ph]]. 

L. Basso, O. Fischer and J. J. van der Bij, Europhys. Lett. 101 (2013) 51004 
[arXiv:1212.5560]. O. Fischer and J. J. van der Bij, JCAP 1401 (2014) 01, 032 
[arXiv:1311.1077 [hep-ph]]. 

[19] P. Fileviez Perez and S. Spinner, Phys. Lett. B 723 (2013) 371 [arXiv:1209.5769 [hep-ph]]. 


36 - 



[20] K. Agashe, A. Azatov, A. Katz and D. Kim, Phys. Rev. D 84 (2011) 115024 
[arXiv:1109.2842 [hep-ph]]. 

[21] CMS-HIG-13-004, arXiv:1408.3316 

[22] G. Aad et al. [ATLAS Collaboration], JHEP 1503 (2015) 088 [arXiv:1412.6663 [lrep-ex]]. 
CMS Collaboration [CMS Collaboration], CMS-PAS-HIG-14-020. CMS Collaboration [CMS 
Collaboration], CMS-PAS-HIG-13-026. 

[23] P. Bandyopadhyay and E. J. Chun, JHEP 1505 (2015) 045 [arXiv: 1412.7312 [hep-ph]]. 

P. Bandyopadhyay, E. J. Chun and J. C. Park, JHEP 1106 (2011) 129 [arXiv: 1105.1652 
[hep-ph]]. P. Bandyopadhyay, E. J. Chun, H. Okada and J. C. Park, JHEP 1301 (2013) 079 
[arXiv:1209.4803 [hep-ph]]. C. Coriano, L. Delle Rose and C. Marzo, Phys. Lett. B 738 
(2014) 13 [arXiv: 1407.8539 [hep-ph]]. 

[24] G. Aad et al. [ATLAS Collaboration], Phys. Rev. D 90 (2014) 5, 052005 [arXiv: 1405.4123 
[hep-ex]]. K. A. Olive et al. [Particle Data Group Collaboration], Chin. Phys. C 38 (2014) 
090001. 

[25] C. Arina, E. Del Nobile and P. Panci, Phys. Rev. Lett. 114 (2015) 011301 [arXiv: 1406.5542 
[hep-ph]]. A. V. Belikov, J. F. Gunion, D. Hooper and T. M. P. Tait, Phys. Lett. B 705 
(2011) 82 [arXiv: 1009.0549 [hep-ph]]. 

[26] F. Staub, Comput. Phys. Commun. 184 (2013) pp. 1792 [Comput. Phys. Commun. 184 
(2013) 1792] [arXiv: 1207.0906 [hep-ph]]. 

[27] A. Pukhov, “CalcHEP 3.2: MSSM, structure functions, event generation, batchs, and 
generation of matrix elements for other packages”, [arXiv:hep-ph/0412191]. 

[28] A. Delgado, G. Nardini and M. Quiros, arXiv:1303.0800 [lrep-ph]. 

[29] P. Franzini, D. Son, P. M. Tuts, S. Youssef, T. Zhao, J. Lee-Franzini, J. Horstkotte and 
C. Klopfenstein et al., Phys. Rev. D 35 (1987) 2883. 

J. S. Lee and S. Scopel, Phys. Rev. D 75 (2007) 075001 [hep-ph/0701221 [HEP-PH]]. 

[30] S. Chatrchyan et al. [CMS Collaboration], Phys. Rev. Lett. 109 (2012) 121801 
[arXiv:1206.6326 [hep-ex]]. 

[31] M. J. Dolan, C. McCabe, F. Kahlhoefer and K. Schmidt-Hoberg, JHEP 1503 (2015) 171 
[arXiv:1412.5174 [hep-ph]]. S. Andreas, O. Lebedev, S. Ramos-Sanchez and A. Ringwald, 
JHEP 1008 (2010) 003 [arXiv: 1005.3978 [hep-ph]]. 

[32] Y. Mambrini, J. Phys. Conf. Ser. 375 (2012) 012045 [arXiv:1112.0011 [hep-ph]]. 

[33] R. Bernabei, P. Belli, F. Montecchia, F. Nozzoli, F. Cappella, A. Incicchitti, D. Prosperi and 

R. Cerulli et al., Int. J. Mod. Phys. A 21 (2006) 1445 [astro-ph/0511262], 

[34] C. Coriano, N. Irges and E. Kiritsis, Nucl. Phys. B 746 77 (2006), C. Coriano, N. Irges and 

S. Morelli Nucl. Phys. B789 (2008) 133. 

[35] C. Coriano, M. Guzzi, and G. Lazarides, Phys. Rev D 82 (2010) 065013. C. Coriano and M. 
Guzzi, Nucl.Phys. B 826 (2010) 87. 

[36] C. Coriano, M. Guzzi and A. Mariano, Phys. Rev. D85 (2012) 095008. C. Coriano, M. 
Guzzi, N. Irges and A. Mariano, Phys.Lett. B 671 (2009) 87. 

[37] P. Bandyopadhyay, C. Coriano and A. Costantini, in preparation. 

[38] V. Khachatryan et al. [CMS Collaboration], Eur. Phys. J. C 74 (2014) 10, 3076 
[arXiv: 1407.0558 [hep-ex]]. 


37- 



[39] G. Aad et al. [ATLAS Collaboration], Phys. Rev. D 90 (2014) 11, 112015 [arXiv: 1408.7084 
[hep-ex]]. 

[40] P. Athron, M. Muhlleitner, R. Nevzorov and A. G. Williams, JHEP 1501 (2015) 153 
[arXiv:1410.6288 [hep-ph]]. 

[41] D. Curtin, R. Essig, S. Gori, P. Jaiswal, A. Katz, T. Liu, Z. Liu and D. McKeen et al ., Phys. 
Rev. D 90 (2014) 7, 075004 [arXiv:1312.4992 [hep-ph]]. 

[42] A. Djouadi, Phys. Rept. 457 (2008) 1 [hep-ph/0503172]. 

[43] A. Djouadi, Phys. Rept. 459 (2008) 1 [hep-ph/0503173]. 

[44] ATLAS-CONF-2013-017 

[45] P. Bandyopadhyay, JHEP 1108 (2011) 016 [arXiv:1008.3339 [hep-ph]]. 

[46] P. Bandyopadhyay and K. Huitu, JHEP 1311 (2013) 058 [arXiv: 1106.5108 [hep-ph]]. 

[47] P. Bandyopadhyay, P. Ghosh and S. Roy, Phys. Rev. D 84 (2011) 115022 [arXiv: 1012.5762 
[hep-ph]]. 



